A direct link between neutrinoless double beta decay and leptogenesis 
in a seesaw model with S4 symmetry 
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We study how leptogenesis can be implemented in a seesaw model with Si flavor symmetry, which 
leads to the neutrino tri-bimaximal mixing matrix and degenerate right-handed (RH) neutrino 
spectrum. Introducing a tiny soft S4 symmetry breaking term in the RH neutrino mass matrix, we 
show that the flavored resonant leptogenesis can be successfully realized, which can lower the seesaw 
scale much so as to make it possible to probe in colliders. Even though such a tiny soft breaking 
term is essential for leptogenesis, it does not significantly affect the low energy observables. We also 
investigate how the effective light neutrino mass |(m ee )| associated with neutrinoless double beta 
decay can be predicted along with the neutrino mass hierarchies by imposing experimental data of 
low-energy observables. We find a direct link between leptogenesis and neutrinoless double beta 
decay characterized by |(m ee )| through a high energy CP phase <j>, which is correlated with low 
energy Majorana CP phases. It is shown that our predictions of |{m ee }| for some fixed parameters 
of high energy physics can be constrained by the current observation of baryon asymmetry. 

PACS numbers: 



I. INTRODUCTION 

Recent experiments of the neutrino oscillation go into 
a new phase of precise determination of mixing angles 
and mass squared differences indicating that the 
tri-bimaximal (TBM) mixing for three flavors can be 
regarded as so-called PMNS mixing matrix C/pmns = 
UtbPv m the lepton sector [2| 
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and P v is a diagonal matrix of phase for Majorana neu- 
trinos. However, properties related to the leptonic CP 
violation are completely unknown yet. The large mixing 
angles, which may be suggestive of a flavor symmetry, are 
completely different from the quark mixing ones. In last 
few years there have been lots of efforts in searching for 
models which produce the TBM pattern for the neutrino 
mixing matrix, and a fascinating way seems to be the use 
of some discrete non-Abelian flavor groups added to the 
gauge groups of the Standard Model. There is a series of 
models based on the symmetry group A4 T' [|| and 
more recently S4 0-0] • 
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In addition to the explanation for the smallness of 
neutrino masses, type-I seesaw model [8|, in which 
heavy right-handed singlet Majorana neutrinos are in- 
troduced, has another appealing feature so-called lep- 
togenesis mechanism for the generation of the observed 
baryon asymmetry of the Universe (BAU) through the 
decay of heavy Majorana neutrinos |9[ . If this BAU orig- 
inated from leptogenesis, then CP symmetry in the lep- 
tonic sector must be broken. So any observation of the 
leptonic CP violation, or demonstrating that CP is not a 
good symmetry of the leptons, can strengthen our belief 
in leptogenesis. For Majorana neutrinos there are two 
additional phases in C/pmnSi one ( or a combination) of 
which in principle can be explored through neutrinoless 
double beta (Oz/2/3) decay [Tg. Although the exact TBM 
mixing pattern forbids low energy CP violation measur- 
able in neutrino oscillation due to U e 3 = 0, it may allow 
CP violation due to the Majorana phases. Therefore, it 
is interesting to explore the existence of CP violation in 
the lepton sector due to the Majorana CP phases in the 
light of leptogenesis. Since |(m ee )| depends on the Majo- 
rana CP phases, we examine if there exists a link between 
0v2f3 decay and BAU. 

In the neutrino models with some discrete flavor sym- 
metries, it is worthwhile to examine if leptogenesis can 
work out while keeping TBM pattern for neutrino mixing 
matrix. Motivated by this issue, in this letter, we study 
how leptogenesis can work in a seesaw model with S4 fla- 
vor symmetry ll]. As anticipated, leptogenesis can not 
work in seesaw models with exact SU(2)l x U(1)y X S4 
symmetry mainly due to the fact that the discrete sym- 
metry S4 leads to zero value of U e 3 and thus could not 
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generate lepton asymmetry. To make leptogenesis suc- 
cessfully realized, it is essential to break S4 symmetry. 
In the case that the right-handed (RH) neutrinos are hi- 
erarchical in mass, successful leptogenesis requires that 
the RH neutrinos are superheavy, which makes inacces- 
sible to colliders and gives rise to the overproduction of 
gravitinos during reheating in the supersymmetric sce- 
narios [12[ . These problems can be avoided if we consider 
an almost degenerate heavy RH neutrino mass spectrum 
so that lepton asymmetry can be resonantly generated 
[l3l IhH ]. In this paper, we especially study how this so- 
called resonant leptogenesis can be implemented in a see- 
saw model with 64 flavor symmetry |15f In this respect, 
the model has to produce degenerate RH neutrino mass 
spectrum at the leading order. The S4 seesaw model rele- 
vant to our purpose has been proposed in [7[ ■ We examine 
how resonant leptogenesis can be implemented by intro- 
ducing a small perturbation in heavy RH neutrino mass 
matrix Mr while keeping Dirac neutrino mass matrix and 
charged lepton mass matrix unchanged in the model Q ■ 
We also show that leptogenesis can be linked to the neu- 
trinoless double beta decay through seesaw mechanism 
with a small perturbation in Mr. 

This work is organized as follows. In Sec. II, we study 
how low energy neutrino oscillation observables are pre- 
dicted in a supersymmetric seesaw model based on the 
flavor symmetry group S4. We also discuss about the 
effective neutrino mass associated with neutrinoless dou- 
ble beta decay. In Sec. Ill, we examine how successful 
leptogenesis can be implemented by introducing a soft 
symmetry breaking term in the heavy RH neutrino mass 
matrix and show how leptogenesis can be liked to 0i>2/3 
decay. Sec. IV is devoted to our conclusion. 

II. LOW ENERGY OBSERVABLES 

Although there have been several proposals to con- 
struct lepton mass matrices in the framework of seesaw 



incorporating S4 symmetry [l6| . in this paper, we 
consider the model proposed in [7[, which gives rise to 
TBM mixing pattern of the lepton mixing matrix and 
leads to degenerate heavy RH neutrino mass spectrum 
which is essential for resonant leptogenesis. The model 
is supersymmetric and based on the flavor discrete group 
Gf = S4 x Z3X Z4, where the three factors play different 
roles. The reasons, why we consider superymmetry, are 
to simplify the choice of desirable vacuum alignment, and 
low scale leptogenesis is well motivated in supersymmc- 
try thanks to the gravitino problem. In general, there ex- 
ists a contribution to leptogenesis via scalar right-handed 
sneutrinos, but their effect depends on soft susy breaking 
terms. In our study, we do not consider the contribu- 
tion by simply assuming that the values of the soft susy 
breaking parameters do not give rise to the contribution 
to leptogenesis. 

The S4 component controls the mixing angles, the aux- 
iliary Z3 symmetry guarantees the misalignment in flavor 
space between the neutrino and the charged lepton mass 
eigenstates, and the Z4 component is crucial to elimi- 
nating the unwanted couplings and reproducing the ob- 
served mass hierarchy. In this framework the mass hi- 
erarchies are controlled by the spontaneous breakdown 
of the flavor symmetry instead of the Froggatt-Nielsen 
mechanism [uj. S4 is the discrete group given by the 
permutations of four objects and has been studied in lit- 
erature but with different aims and different results. 
It is composed by 24 elements, divided into 5 irreducible 
representations: two singlets, Ii and I2, one doublet, 
2, and two triplets, 3i and 32. The technical details of 
the group are shown in Q. The multiplication rules be- 
tween the various representations are as follows; li©Ij = 

mod 2) + l; lj <8> 2 = 2, lj ® 3j = 3((j+j) mo d 2)+l> 

2®2 = liel 2 ©2,2(g)3 i = 3i®3 2 , 3*® 3* = li©2©3i©3 2 , 
3i ® 3 2 = 1 2 © 2 © 3 2 © 3 2 , with i,j = 1, 2. The matter 
fields in the lepton sector and the flavons under Gf of 
the model are assigned as Table Q] 



TABLE I: Representations of the matter fields in the lepton sector and the flavons under S4 x Z3 x Z4. And ui is the third root 
of unity, i.e. u = e i2 * /3 . 
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The superpotential of the model in the lepton sector reads as follows 
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where the dots denote higher order contributions. 

The vacuum configuration can be determined from the 
vanishing of the derivatives of the superpotentials wi and 
W v with respect to each component of the driving fields, 
as shown in [3] . Using this way, we can obtain the align- 
ment of the VEVs of flavons as follows; 

{<p) = (0, v v , 0) , ( X ) = (0, v x , 0) , (0) = w , 
(v) = (v v ,v v ) , (cf>) = (vfaVfaVf) , (A) =v A .(4) 

With these VEVs alignments as well as breaking of elec- 
troweak gauge symmetry, the charged-lepton mass ma- 
trix is explicitly expressed as 
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where we assume all components are real, and the neu- 
trino sector gives rise to the following Dirac and Majo- 
rana matrices 
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FIG. 1: Allowed parameter region of the ratio r — b/a as a 
function of cos <j> constrained by the la experimental data in 
Eq. (|13[) . The thickness of the line reflects la uncertainty due 
to experimental data. Here, the blue (dark) and red (light) 
curves correspond to the inverted and normal mass ordering 
of light neutrino, respectively. 



with r = b/a, and too = ^a 2 /M where v u = usin/3 
and v = 176 GeV. And the lepton mixing at low energy 
t/pMNS = U v is given by 



/ 2be l * a - be** a - be 1 * 



C^PMNS 



U TB V>i&g.{l,e 1 ^ 



JP2 



(10) 



Mr 



a -be** a + 2be** 
\a- be** -be 1 * 

[moo 

M 
y M 



-be** 
a + 2be** 



(6) 



(7) 



where we assume the quantity M, a and b are real and 
positive quantities, and the relative phase <fr = a 2 — a% 
is the only physical phase because the phase a\ can 
be rotated away. After seesawing, the effective light 
neutrino mass matrix is obtained from seesaw formula 
TO c ff = — {m^) T M^m^, which can be diagonalized by 
the TBM mixing matrix as follows; 



Ujm ef sU v = Diag.(TOi,m 2 ,TO 3 ), 

where the mass eigenvalues are given as 

toi = too(1 + 9r 2 — 6r cos</>) , 
to 2 = 4too , 

«3 = toq(1 + 9r 2 + 6rcos</>) , 



(8) 



(9) 



where (i\ — -fi/2, (3 2 — (71 — 72)/2 are the Majorana 
CP-violating phases, with 

71 = arg{(a-3&e^) 2 } , 

72 = arg{-(a + 36e^) 2 } . (11) 

The phase factored out to the left have no physical mean- 
ing, since it can be eliminated by a redefinition of the 
charged lepton fields. There are the nine physical quan- 
tities consisting of the three light neutrino masses, the 
three mixing angles and the three CP-violating phases. 
The mixing angles are entirely fixed by the Gf symme- 
try group, predicting TBM and in turn no Dirac CP- 
violating phase, and the remaining 5 physical quantities 
Pi, /3 2 , toi, m 2 and to 3 , are determined in terms of the 
five real parameters M, v u ,a,b and </>. 

Because of the requirement of MSW resonance for so- 
lar neutrinos implying Ato 2 ^ > 0, and the observed hi- 
erarchy IAtoJ-jJ 3> Am^x, there are two possible neu- 
trino mass ordering depending on the sign of cos </>: (i) 
m i < m 2 < m 3 (normal mass ordering) which corre- 
sponds to cos (f> > and (ii) to 3 < mi < m 2 (inverted 
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FIG. 2: Predictions of the effective mass |(m ee )| for 0^2/3 as a function of the phase <f> in the left panel and the ratio r in the 
right panel based on the la experimental results given in Eq. (|13[) . Here, in both panels the red (light) and blue (dark) curves 
correspond to the normal mass spectrum of light neutrino and the inverted one, respectively. 



mass ordering) which corresponds to cos (f> < 0. The so- 
lar and atmospheric mass-squared differences, which are 
given by 

Am^ = 3mo(l - 3r 2 -I- 2r cos </>) (5 + 9r 2 - 6rcos0) , 
\Am 2 31 \ = 24m 2 r|cos< ? (>|(l + 9r 2 ) , (12) 

are constrained by the neutrino oscillation experiments. 
Since neutrino oscillation data indicate that Am^ is pos- 
itive, 1 — 3r 2 + 2r cos <fi > 0. It is interesting to see how 
the ratio Am^/ Arn^j leads to the correlation between 
the paremeters r and cos(/>. From Eq. (fl"2"j) . we see that 
the ratio Am^/Am^ is independent of mo, and that 
cos0 = is not allowed, which is reflected in Fig. 1. For 
our purpose, we consider the experimental data at 1 a [1|: 

| Am 2 .! | = (2.29 - 2.52) x !CT 3 eV 2 , 



Am 2 .! = (7.45 - 7.88) x l(T 5 cV 2 . (13) 



Hereafter, we use the la confidence level experimental 
values of low energy observables for our numerical cal- 
culations. Imposing the experimental results, we present 
the correlations between r and cos cj> for normal mass hi- 
erarchy (red-plot) and inverted one (blue plot) in Fig. [1] 
From Eqs. (|12I13I) . we can also determine the value of 
mo as a function of cos 4>. Thanks to the zero entry 
in t/pMNS, m 3 does not contribute to the effective neu- 
trino mass and thus only the phase /3\ can contribute 
to the 0v2(3 decay amplitude, which can be written as 
|(m ee )| = \mi\U e i\ 2 + m 2 \U e2 \ 2 + m 3 \U e3 \ 2 \, where U el 
(i=l,2,3) are the components of C/pmns- As a result, the 
effective mass governing the 0v2(3 decay is given by Q 



|(m ee )| = i|2m ie 2 ^+m 2 | 



2m y/l -4rcos0 + 2r 2 (2 + 3cos2</>) - 12r 3 cos </> + 9r 4 
I 



(14) 



The behavior of |(m ee )| is plotted in Fig. [2] as a func- 
tion of the phase <fi (left panel) and the ratio r (right 
panel). Here we note that mo is related with other pa- 
rameters through the seesaw formula mo = a 2 v 2 /M. For 
our numerical analysis, we take mo = 0.0035—0.04, mo — 
0.012 — 0.05 for normal and inverted orderings, respec- 
tively, which are consistent with the experimental results 
given in Eq. (1131) . Numerically, our prediction is turned 
out to be 0.0076eV < |(m ee )| < 0.10 eV for the normal 
mass ordering (red-plot) and 0.043eV < |(m ee )| < 0.12 
eV for the inverted mass ordering (blue-plot), where the 
upper limits come from the cosmological bound on neu- 
trino mass scale. We note that the existence of the lower 
bound on the prediction of m ee in the case of the nor- 



mal mass ordering is due to |cos<jf>| < 1, which is easily 
understood from Eq. ©, as also discussed in Q- In 
Fig. [2 the horizontal solid and dashed lines correspond 
to the current lower bound sensitivity (0.2 eV) Jl8j and 
the future lower bound sensitivity (10 -2 eV) [19| of 0^2/3 
experiments, respectively. The thickness of the lines cor- 
respond to the Id allowed ranges due to experimental 
results. 

Using Eq. ([TTj) we can obtain the explicit correlation 
between the phase <p and the Majorana phase fi\ 
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Note here that the size of tan2/?i is constrained by Fig. 
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FIG. 3: Correlation of the Majorana CP phase /3i with the 
phase <j> constrained by the la experimental data in Eq. (|13[1 . 
The red (light) and blue (dark) curves correspond to the nor- 
mal mass spectrum of light neutrino and the inverted one, 
respectively. 

[3] represents the correlation between the phase <f> and the 
Majorana phase /3i for normal mass ordering (red-plot) 
and inverted one (blue-plot). 

III. LEPTOGENESIS THROUGH SOFT S A 
BREAKING 

In a basis where the charged current is flavor diagonal, 
the diagonalization of Mr leads to the mass eigenvalues 



of heavy RH neutrinos given by 

V^M R V R = Diag.(M, M, —M) , (16) 
where the diagonalizing matrix Vr is 

v «=\ 73 7f ]• 

After performing a basis rotation so that the heavy RH 
Majorana mass matrix Mr becomes diagonal by the uni- 
tary matrix Vr, the Dirac mass matrix mf, gets modified 
to 

mt -> Y v v u = V%m d u (17) 
where the Yukawa coupling matrix Y v , is given as 

(26e^ a - be^ a - be i(t > \ 

V2(a - be 1 *) ^ff- 2±^1 _ (18) 

n q+3be'^ q+36e^ / 

U V2 V2 I 

Here, we notice that the CP phase <j> existing in mf, obvi- 
ously gives rise to low-energy CP violation. On the other 
hand, leptogenesis is associated with both Y v itself and 
the combination of Yukawa coupling matrix, H = Y U YJ , 
which is given as 



H = 



I 2a 2 + 6b 2 - 4ab cos 4> 
V2{a 2 - 36 2 + 2a6cos< 

V 



V2{a 2 - 36 2 + 2a6cos< 
3a 2 + 3b 2 - 2abcos(j) 







% 2 + 6ab cos < 



(19) 



We see that all Im[7Jy ■] vanish and thus leptogenesis could 
not take place. However, when we consider low scale 
much below 10 9 GeV to avoid the abundance of gravitino 
overproduction problem [22j within a super-symmetric 
version, the r and /x Yukawa couplings comes into equi- 
librium as the lepton asymmetry is being created. In 
that case, lepton number asymmetries and washout ef- 
fects become flavor dependent, and this can give rise to 
a final baryon asymmetry which is different in size from 
the one generated if flavor effects are ignored [13, Hl[ . 
We call it flavored leptogenesis that occurs in general for 
T ~ Mj(l + tan 2 0) < 10 9 GeV in supersymmetric see- 
saw model [H] ■ In addition, since exact degenerate mass 
spectrum of the heavy RH neutrinos prevents leptogene- 
sis to be occurred, in order for flavored leptogenesis to be 
viable, the degeneracy of the heavy RH neutrino masses 
should be lifted. 



For an almost degenerate heavy Majorana neutrino 
mass spectrum, leptogenesis can be naturally imple- 
mented through the resonant-leptogenesis framework 
Guam. In this case, the CP asymmetry generated by 
the i-th heavy Majorana neutrino decaying into a lepton 
flavor a is given by [241 ] 



Im \h vj (Y^ia (Y„)* Q 







*1 



AM^ 



(20) 



where Tj = HjjMj/8ir is the decay width of the j-th 
right-handed Majorana neutrino and 5% is the mass split- 
ting parameter defined as 



N 



1 



Mi 



(21) 



In order for resonant leptogenesis to be implemented 
successfully, soft terms of the form eM vfv? with small 
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dimensionless parameter e [25j are added in Eq. (3), 
which lift the degeneracy of the heavy Majorana neu- 
trino masses. Although there are several possibilities to 
incorporate the breaking parameter e in Mr, which lead 
to the mass slitting parameter \6^\ ~ e. Without loss 
of generality, we introduce a breaking term of the form 
tMv^y^, which modifies the RH neutrino mass matrix as 



f M 
eM M 
y M 



(22) 



where the parameter e is assumed to be real. Mr is di- 
agonalized as V^MrVr = Diag.(Mi, M 2 , M 3 ) with real 
eigenvalues given as 



Mx = M , M 2 ^M(1+-} , M 3 



-M(1-~J (23) 



and the diagonalizing matrix Vr is written as 
/ 1 



Vr 



n J- + ^1 -I f— 



(24) 



y/2 4\/2 4%/2 



+ 



In the basis where the heavy RH mass matrix real and 
diagonal, the Dirac Yukawa coupling matrix now reads 
as 



1 ~ T 
Y» = —Vr? 



(25) 



where mf, and Vr are given in Eq. ^ and Eq. 
respectively. Assuming e C 1, we find that the ma- 
trix Y V YJ is real and almost the same as Y^YJ ~ H 
given in Eq. (TH)1) . As a result, the contributions of 
N3 to lepton asymmetries ef can be negligible, due to 



Hi 
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3(31) 



H. 



23(32) 



jorana neutrino masses given in Eq, 
the mass splitting parameter as follows, 
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N 



°N 



0. From the heavy Ma- 
we can obtain 



(26) 



Then, combining with Eqs. d [19]) and Eq. {26]), the 
flavor dependent CP asymmetries ef can be obtained as 
follows 



a 2 (I + 2rcos0 - 3r 2 



e47r(l — 2r cos ( 
a 2 (l + 2r cos</ 



-3r 2 ) 
3r 2 ) 



e27r(3-2rcos0 + 3r 2 ) 



a 2 (l + 2rcos0-3r 2 ) 
£8tt(1 - 2rcos<?!) + 3r 2 ) ) 
a 2 (l + 2rcos0-3r 2 ) 
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e47r(3 — 2r cos + 3r 2 ) 



(27) 



In these expressions, the values of the parameters <j> and 
r are determined from the analysis as demonstrated in 
sec. II, whereas e and a are arbitrary. However, since the 
seesaw relation a 2 — m M/v 2 as defined in Eq. ©, the 
value of a depends on the magnitude of M once m is 
determined. Thus, in our numerical analysis, we take M 
and e as input in the estimation of lepton asymmetry. 
Here, we note that although e and M are taken to be 
independent parameters in our analysis, the predictions 
of the lepton asymmetries e" 2 depends only on the ratio 
M/e. We see from Eq. ([27]) that the asymmetries can be 
substantially enhanced by lowering e. However, the value 
of e is constrained by the validity of the perturbation. In 
order for the perturbative approach to be valid, the tree- 
level decay width I\ must be much smaller than the mass 
difference: 

T 4 = ^^" Aij < M 2 - M x = S N M 2 , i = 1, 2 . (28) 

SIT 

It is approximately given by 

r, ; ~ a 2 (l + 3r 2 - 2rcos</>)/47r < S 2 }. 

Due to a 2 = m M/v 2 and d 2 } ~ e/2, if we take the 
seesaw scale, as an example, to be M = 10 6 GeV, the 



value of a is the order of O(10 5 ), which requires e 3> 
1(T 10 for (5at ^ e/2. 

Once the initial values of ef are fixed, the final result 
of rjB can be obtained by solving a set of flavor dependent 
Boltzmann equations including the decay, inverse decay, 
and scattering processes as well as the nonperturbative 
sphaleron interaction. In order to estimate the wash-out 
effects, we introduce the parameters Kf which are the 
wash-out factors due to the inverse decay of the Majorana 
neutrino iVj into the lepton flavor a(— e, fj,, r) [26]. The 
explicit form of Kf is given by 

pa ^2 

= t_7Y lur \ = (Yv )oti(Yi/)ia~ TT (29) 
H(Mi) m^Mi 

where Tf is the partial decay width of Ni into lepton fla- 
vor a and Higgs scalars, H(Mi) ~ {^g*/45)?M 2 jM Pl 
with the Planck mass M Pl = 1.22 x 10 19 GeV and the 
effective number of degrees of freedom ~ 228.75 is the 
Hubble parameter at temperature T = Mi , and the equi- 
librium neutrino mass ~ 10~ 3 . From Eqs. (|18|) and 
([29]) we can obtain the washout parameters as follows 

K{ * 4r 2 ^ , 
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FIG. 4: The prediction of r\B as a function of |(m ee )| for M = 10 6 GeV and the soft breaking parameter e = 10 -7 . The left-plot 
and the right-plot represent the normal mass ordering of light neutrino and the inverted one, respectively. The solid horizontal 
line and the dotted horizontal lines corresponds to the experimental value of baryon asymmetry, rjg MB = 6.1 x 10 -10 , and 
phenomenologically allowed regions 2 x 10~ 10 < i]b < 10 -9 , respectively. 



mo , ON 
— (1 -2rcos0 + r 2 ) , 

TO* 

^^(1 - 2rcos0 + r 2 ) 

TO* 

-^-(1 + 2r cos (j) + r 2 ) 
2to* 



(30) 



Since we take M to be 10 6 GeV, each lepton asymmetry 
for a single flavor in Eq. (j2"7| is weighted differently by 
the corresponding washout parameter given by Eq. (|30[) , 
and appears with different weig ht in the final formula for 
the baryon asymmetry pU. |26| as follows; 
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r 



with wash-out factor 

„ /8.25 



KT 



V0.2 



1.16 



(32) 



As can be seen from Eqs. (|27|30j) . since the lepton asym- 
metries in fi and r flavors are equal to the first order, 
satisfying £^ 2 , +e[( 2 ) = ~ £ i(2)j an{ i the washout factors 
in \x and r are also equal, the value of baryon asymmetry 
can be obtained as 



i lB ~ 10- 2 {e?« - K\) + e%{k» - «|)} 



(33) 



Numerically, the washout factors for normal mass spec- 
trum of neutrino and the inverted one are given as 



< 0.06, 



< 0.05 



< 0.02, «£-«§< 0.05 , 



(34) 



respectively. Thus, taking Eqs. (|33I34[) into account, 
is needed to obtain a successful leptoge- 



-1(2) 



10 
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nesis, which in turn means that the value of soft breaking 
parameter e < 10 -6 is required for M = 10 6 GeV. 

Before going to discuss the value of tjb, we consider 
the effects of soft breaking term on the light neutrino ob- 
servables. The effective neutrino mass matrix is modified 



due to the soft breaking term and thus it can be diago- 
nalized as U^m c eU u = Diag.(mi, TO2, TO3) with the real 
eigenvalues 



fhi = rtii + m,0(e) , 



(35) 



and the mixing matrix U v = C/pmns can be written as 



PMNS 



UtbPu + 5U{e)P v 



(36) 



) is diagonal matrix of Ma- 



where P v — Diag.(l, 
jorana phases defined in Eq. (fTUj) . Note here that, since 
e is real, the Majorana CP phases (3\^ are not affected 
from the soft breaking. However, it is obvious that the 
effects of e to the light neutrino mass eigenvalues and the 
neutrino mixing angles are negligible due to e < 10~ 6 . 
As expected, the value of |(m ee )| is almost the same as 
Eq. (PI]). 

The predictions for r]B as a function of |(m ee )| are 
shown in Fig. @] where we have used M = 10 6 GeV, 
e = 10~ 7 and tan/3 = 2.5 as inputs. Please note that 
if there is a mass splitting of Majorana neutrinos by a 
renormalization group equation, the value of tan /? can 
be crucial to have a successful leptogenesis. However, 
its value does not affect significantly our results. So, for 
simplicity, we take tan/3 = 2.5. The horizontal solid and 
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dashed lines correspond to the central value of the experi- 
ment result of BAU ryg MB = 6.1 x 10" 10 H3 and the phe- 
nomenologically allowed regions 2 x 10 -10 < t\b < 10~ 9 , 
respectively. As shown in Fig. 2J the current observa- 
tion of f7§ MB can narrowly constrain the value of |(m ee )| 
for the normal hierarchical mass spectrum of light neu- 
trino and inverted one, respectively in the case that the 
scale of leptogenesis is 10 3 < M < 10 6 GeV and e is small 
enough. Combining the results presented in Figs. [2] and [3] 
with those from leptogenesis, we can pin down the Ma- 
jorana CP phase /3i via the parameter <f>. 

IV. CONCLUSIONS 

We have examined how leptogenesis can be imple- 
mented in the seesaw model with S4 flavor symmetry 
which leads to the neutrino TBM mixing matrix and de- 
generate RH neutrino spectrum. In order for leptogenesis 
to be viable, Y v Yj should contain nontrivial imaginary 
part and the degeneracy of heavy right-handed Majorana 
masses has to be lifted. In our study, we have shown that 
flavored resonant leptogenesis can be successfully realized 
by introducing a tiny soft S4 symmetry breaking term 
in Mr, which can much lower the seesaw scale. Even 
though we can lower the seesaw scale down to TeV scale, 
it would be difficult to probe it directly in collider exper- 
iments because the couplings between the heavy neutri- 
nos and the light neutrinos are so small in our scenario. 
Even though such a tiny soft breaking term is essential 



for leptogenesis, it does not significantly affect the low 
energy observables. We have also investigated how the 
effective neutrino mass |(m ee )| associated with 0^2/3 de- 
cay can be predicted along with the light neutrino mass 
hierarchies in our scenario for leptogenesis by imposing 
experimental data of low-energy observables. Interest- 
ingly enough, we have found a direct link between lepto- 
genesis and neutrinoless double beta decay characterized 
by |(m ee )| through a high energy CP phase <j> which is 
correlated with low energy Majorana CP phases. We 
also have shown that our predictions of |(m ee )| for some 
fixed parameters can be constrained by the current ob- 
servation of baryon asymmetry in the case that the scale 
of leptogenesis is 10 3 < M < 10 6 GeV and e is small 
enough. 
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